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ABSTRACT 

Recent work has shown that the speed of the cooling front in soft X-ray 
transients may be an important clue in understanding the nature of accretion 
disk viscosity. In a previous paper (Vishniac and Wheeler 1996) we derived 
the scaling law for the cooling front speed. Here we derive a similarity solution 
for the hot inner part of disks undergoing cooling. This solution is exact in 
the limit of a thin disk, power law opacities, and a minimum hot state column 
density which is an inhnitesimal fraction of the maximum cold state density. 
For a disk of finite thickness the largest error is in the ratio of the mass flow 
across the cooling front to the mass flow at small radii. Comparison to the 
numerical simulations of Cannizzo et ah (1995) indicates that the errors in other 
parameters do not exceed {csp/rF^Fy, that is, the ratio of the sound speed at 
the disk midplane to its orbital velocity, evaluated at the cooling front, to the 
gth power. Here g ~ 1/2. Its precise value is determined by the relevant hot 
state opacity law and the functional form of the dimensionless viscosity. 


1. Introduction 

The most popular model for soft X-ray transients and dwarf novae is that both are due 
to fluctuations in the luminosity from accretion disks surrounding black holes and white 
dwarfs, respectively, and that the specific mechanism that starts an outburst is a thermal 
instability associated with the ionization of hydrogen (for a recent review see Cannizzo 
1993a or Osaki 1996). In addition to their intrinsic interest, these systems may provide 
valuable clues to the nature of angular momentum transport in disks, simply by virtue 
of being non-stationary systems. In the thermal instability model these disks make the 
transition back to their quiescent states as their outer parts cool and a thermal transition 
front propagates inward. Mineshige, Yamasaki, and Ishizaka (1993) pointed out that this 
model will produce the observed exponential decay of soft X-ray luminosity for soft X-ray 
transients only if the radius of the hot phase decreases exponentially, i.e. if the cooling front 
speed is proportional to the radius of the uncooled part of the disk. 
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The speed of the cooling front affects the evolution of the hot phase of the disk and its 
value will depend on the local physics of the disk, including the dimensionless viscosity a. 
More specihcally, Mineshige (1987) and Cannizzo, Shatter and Wheeler (1988) argued that 
the cooling front speed is approximately 


Tp ~ -OLpCsF 




( 1 ) 


where ap is the dimensionless viscosity at the onset of rapid cooling, H is the disk thickness, 
Sr is the width of the cooling region, c* is the sound speed in the midplane of the disk, 
and the subscripts ‘F’ denote values at the onset of rapid cooling. Numerical simulations 
seemed to show that Sr ~ O.lr. 


Recent high resolution work by Cannizzo, Chen and Livio (1995, hereafter CCL) 
has shown that this estimate for Sr was a consequence of insufficient resolution in the 
simulations. They found Sr was approximately equal to the geometric mean of the disk 
height and radius. Combining this with equation (|l|) they argued that an exponential 
decline implied a scaling law for a of the form 


a ~ 50 


CsF 
r pVlp 


3/2 

OC 



( 2 ) 


More recently Cannizzo (1996) has shown that a similar result can be recovered from 
models of the dwarf nova SS Cygni. 


In a previous paper ([Vishniac fc Wheeler 1996| , hereafter VW) we showed that the 
cooling front speed was actually independent of the structure of the disk outside the radius 
where rapid cooling sets in, and that the cooling front speed scales as 


rp OC —OipCsP 



( 3 ) 


where q is determined by the functional form of a and the opacity law for the hot part 
of the disk. Assuming that the hot phase is thermally stable implies that q lies in the 
interval [0,1]. For most reasonable opacity laws g ~ 1/2. This implies that the exponential 
luminosity decline of dwarf novae and soft X-ray transients is not just another test of the 
form of a, but one that specihcally depends only on the physics of the hot state, thereby 
eliminating uncertainties regarding opacities in the cold state and the effects of comparing 
a largely ionized hot state to a largely neutral cold one. Finally, in VW we managed to 
show that the cooling front speed is surprisingly insensitive to the value of a right at the 
cooling front, but instead is a measure of its behavior over a broad annulus inward from the 
cooling front. 
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While it is important to see that the scaling law reported in CCL can be recovered 
from the simple physical argument presented in VW. It would be better still to see that 
the numerical coefficients found by CCL can be reproduced by an analytic solution built 
around this same argument. Furthermore, the numerical solution is not precisely self-similar 
and by itself the scaling law does not give us an understanding of when it might fail. In 
this paper we present a self-similar analytic solution of a cooling wave in a thin accretion 
disk and estimate the error involved in applying it to a disk of hnite thickness. Section 2 
contains a derivation of the thin disk similarity solution. Section 3 contains a summary of 
our results, a discussion of the error involved in applying this solution to realistic disks, and 
a discussion of the implications of this work for the physics of accretion disks. 


2. Analysis 

The usual sign that a system evolving in space and time has a stable self-similar 
solution is that the spatial distribution of some dynamically signihcant variable approaches 
a self-similar shape in an experiment or in a numerical simulation. In the case of a cooling 
wave in a disk the column density in the bulk of the hot phase has a simple power law shape 
with a constant exponent (CCL). However, near the cooling front itself, in the so-called 
‘precursor’ region (cf. VW), the density drops quickly below an extrapolation of this power 
law. In VW we showed that the ratio of the column density at the edge of the precursor 
region, Sp, to the density at the cooling front itself, is given by 

Sp 

Sf \ CsF j 

The fact that this ratio varies with r^, and therefore with time, suffices to show that the 
shape of the column density distribution is not time-invariant and that any self-similar 
solution to this problem can only be an approximation. 

On the other hand, since this column density ratio goes to inhnity as the disk becomes 
inhnitely thin, this result suggests that it may be useful to approximate the column density 
at tf as zero, while retaining a hnite mass outhow determined by the actual conditions at 
the cooling front. We will see that this will allow us to recover an approximate similarity 
solution. This similarity solution is, in a sense, imaginary for the region where S < Sf, but 
we can estimate the errors this procedure introduces by examining the dynamical equations 
near the cooling front. To put it another way, in VW we divided the hot part of the disk 
into an interior region, in which the column density is well-approximated by power laws in 
radius and time, and a narrow precursor region in which the density drops precipitously as 
the cooling front approaches. Here we approximate the entire hot state with a similarity 
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solution, with a thin outer edge which is unrealistic. However, we can use this solution to 
calculate the evolution of the hot state, and use our knowledge of the precursor region to 
understand the errors in this solution. 

One drawback to this approach is that it is not always appropriate to ignore the 
effect of the colder parts of the disk in calculating the cooling front speed and structure. 
More precisely, material ejected across the cooling front will quickly cool and slow down. 
Conservation of mass implies that the column density of cold material just outside the 
cooling front will be comparable to the column density of hot gas at the edge of the 
precursor region. Since there is a maximum column density for the cold state, this implies 
that there is an upper limit to Sp/E^, which is a function of radius and which is comparable 
to l^h,min at that radius. As a cooling wave approaches this limit its speed will drop 

and its structure will become dependent on the details of the cold state. We will ignore this 
limitation in what follows, but it has to be taken into account when the solution presented 
here implies Ep/Ei? very large. We show at the end of this paper that for the systems 
considered so far this does not appear to be an important problem. 

The basic equations for mass flow in a thin disk are 


5tE = -^dr{M), 

2nr 

(5) 


(6) 

M = 27rrEI4-5 

(7) 


where E is the gas column density, Vj. is the radial velocity, Cg is the sound speed, a is the 
dimensionless viscosity, and M is the mass flow. 

The thermal structure of an optically thick disk is determined by its opacity source. 
When the opacity law has a simple power law form, the midplane temperature can be 
written as 

T = (8) 

In the hot portions of the disk, for which the opacity can be approximated by a Kramers 
law opacity, a = c = 3/7 and b = 1/7. Rapid cooling sets in for temperatures below 
where 

( 9 ) 

(cf. CCL, equation (4)). Consequently, Tmin has a very weak dependence on radius. 
Equations and (|^) also imply that the hot phase has a minimum column density 
Emin = Ei?. However, Emin has a strong dependence on radius. 




5 


In what follows we will assume that the dimensionless viscosity a has the form 


a = ao 



( 10 ) 


where n ~ 3/2 is required to get consistency with the nearly exponential decay of flux 
from soft X-ray transients following an outburst. This is not the only form for a that will 
produce an exponential luminosity decay. The luminosity evolution is primarily sensitive to 
the variation of a with radius and the same result can be obtained from a oc ([Vishniac 
&; Wheeler 1996[ ) but equation (|l^) is is the only functional form for a which is consistent 
with all available phenomenological constraints. A more troublesome question is whether 
or not a can be written as a purely local function or whether it depends in some manner 
on the global state of the disk. 

We are looking for an asymptotic state in which the interior region of the disk maintains 
a self-similar mass distribution approximately characterized by a mass accretion rate Mo(t) 
which varies slowly with time. We can write this as 


M{r,t) = Mo{t)F{y), 


( 11 ) 


where 


y = 


rp 


( 12 ) 


At small radii the disk will be very close to an equilibrium state with an approximately 
constant M = Mo{t), so A(0) = 1. At large radii the mass flux will be positive, that is 
outward, towards the cooling front, so F{1) will be negative. Using equations (||) and (pUj) 
we can dehne a constant Cq such that 


aSc? = CnE‘''’r; 


(13) 


where 


and 


1 

Q 


1 + 





1 — c 

T^b- 


(14) 

(15) 


Our choice for the form of equation (|^ is motivated by our previous work in VW in 
which we showed that the mass transfer rate everywhere between the cooling front and the 
inner edge of the disk will vary only by a factor of order unity, although it will change sign. 
This is fairly obvious for the precursor region, since the mass entering this narrow region 
as the cooling front advances has to be approximately balanced by the mass ejected across 











6 


the cooling front. The fact that the asymptotic state of the cooling front is one in which 
the front advances slowly enough to allow the interior region of the disk to evolve, implies 
that the velocity of the cooling front is directly tied to the accretion velocity in the outer 
parts of the interior solution. This in turn implies that the mass transfer rate inward in the 
interior solution is proportional to the mass transfer rate outward through the cooling front. 
It follows that a general description of the cooling disk is most easily obtained by referring 
the mass transfer rate at all radii to its value at the inner edge. (The mass transfer rate 
across the cooling front would be an equally good choice.) 


Using these dehnitions in equation we hnd that 

-rpQpMo fy 


S = 


bvrCn 


-y 


m- 


y 


1/2 


This expression can be substituted into equation to obtain 


gS In Mq + dtlnrp 
It is convenient to rewrite this as 






dy ' 
1/2 


.2 llF{y)r^lHy}} 2nrl 

2y^rp(y) ■ 




y 


1 + C'l 1 


jyF{y)y-^Fdy^ 


= d F 

Uyr , 


where 


and 


dtln rp 
2dt\nMo' 


C 2 = -a* In Mo 


/ 2jrqrl \ 

\ -M„ J 


/'-MoQfrpy 

1 6irC„ J 


(16) 

(17) 

(18) 

(19) 

( 20 ) 


Our assumption that the hot phase of the disk can be described by a self-similar 
solution is equivalent to requiring that C'l and C 2 be constant, so that equation (|^ is 
time-invariant. In fact, Ci can be determined directly if we know the scaling of rp with 
Tp and use the requirement that the various time-varying factors in C 2 combine to give 
a constant product. We will return to this point later. The value of C 2 is determined 
by requiring that the solution to equation (p!8D satisfy the boundary conditions. In other 
words, it is an eigenvalue of equation (p!8[). 


What are the outer boundary conditions for the similarity solution? The condition 


that the column density vanishes at rp implies, using equation ([T^) , that 

dy 


m- 


y 


1/2 


= 0. 


( 21 ) 
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The mass flow at the outer edge of the disk, Mp, is less precisely determined, since it 
depends on the gradient of S at the cooling front. However, given that the cooling front is 
defined by the onset of extremely rapid cooling, it follows that the radial scale for changes 
in E will be close to the disk height and we can write the radial velocity at the cooling front 
as 

Vp = AapCsPi (22) 

where A is a constant of order unity. Here we are interested in the mass flux at an 
imaginary point where E = 0, but we can approximate it with the value at E = Ej?. In 
CCL’s simulation they found A 1/6, but the exact number will depend to some extent on 
the uncertain turbulent radial transport of heat near the cooling front. In any case we have 

Mp = 271T pYjp/^oipCgp. (23) 


This is not a constraint on F{1), which is the ratio between this mass flux and the 
undetermined value of Mo(t). However, we can use this result to rewrite C 2 . Starting with 
the definition in equation (|^) we have 


Co 


We see from this that 


= —dt In Ml 

-a^lnr^^ 


f 2iTqr% 


A 


Mo 


27njr|,Si? 


-Mo 

_ 1 / 2'KqrprF^F 

Mo 


—Mo^ft| 
GttCo 

( —MqQp 
y6napT.pc^^p 


A / —/^rpQp \ ^ 
) {scspFil)) 


2Ci V -. 


rp = 


-6CiC2 / -a y 


\rpyip J 


(24) 


(25) 


in agreement with the scaling derived in VW. 

In order to evaluate Ci from its definition in equation (p!9D we use the fact that in the 
asymptotic state Mp has a fixed ratio to Mq, so that 


Mn = 


27irpAapCspT,p 

F(I) ■ 


Using equation (H) and our assumed form for a we And that 

V rp{l-b{n/2))/a (c-b{n/2))/a 

MK. 1 P Ip , 


(26) 


(27) 


Mq OC yd“^C/2))/g+c-l)/ajqi+fen)/2+(l-fe(n/2))/a 


SO that 


(28) 














Substituting this into the dehnition of Ci we hud that 


^ ^ ^ _ 

2 ((1 - bj) + (c - l)q + €q(a(l + bn)/2 + (1 - bn/2))^ ’ 

where 

Tp OC Tp. 

We can hnd the value of e by using equation (^. It is 

9 - lln 
‘ ^ 140 +lln’ 


( 29 ) 


(30) 


(31) 


which is typically a small number. We note that if n = 3/2 and we can assume a Kramers 
law opacity, then e = —0.048 and Ci = 0.2109. This value of n is consistent with the 
observed approximately exponential luminosity decay, although any value near it would be 
as well. 


If we integrate equation (pK]) for the case with n = 3/2 and Kramers opacities, and 
determine C 2 by requiring that / satisfy the boundary condition given by equation 
then we hnd that 

Mf = - 2 . 88 M 0 , (32) 

f’F = —0.9AaFCsF ^ , (33) 

and M = 0 at r = O.dOrj?. The cooling front velocity goes as Here we have used 

A 1/6 (from CCL). CCL obtained Mp ~ 2Mo and found a zero in M at r Ri 0.38ri?. 
Assuming that the numerical solutions of CCL are more accurate than the similarity 
solution, the errors in the similarity solution are of order 5 to 10% for an example in 
which {H/rY is slightly less than 0.1. The only exception is in the value of F{1), which 
seems to be too large by as much as 40%. For illustration we show S, M, and the surface 
temperature, T^, as a function of r for the similarity solution in hgure 1 . 


3. Error Evaluation and Conclusions 

We have found a similarity solution that describes the progress of a cooling wave in an 
inhnitely thin accretion disk whose hot state is governed by power law opacities. For a disk 
of hnite thickness, the nonzero value of the critical column density for the onset of rapid 
cooling introduces an error at the outer edge of the similarity solution. The free parameters 
of this solution are almost entirely determined by the physics of accretion and the hot state 
opacity, the sole exception being the parameter A, which is determined by the structure 





9 


of the rapid cooling zone. In practice, this means that at present A is determined from 
numerical simulations. We can see, by comparison with the numerical simulations of CCL, 
that if we take a Kramers opacity and n = 3/2 we get a solution which is a fair description 
of the progress of a cooling wave in a disk surrounding a soft X-ray transient. Looking at 
our results in hgure 1, we can see that the temperature prohle of the hot phase is almost 
a perfect power law, with a sharp cutoff near the outer edge. Clearly, it will be difficult 
to observe the cutoff structure directly. However, this model predicts that the power law 
temperature prohle of the inner region will be well above the critical temperature for the 
thermal instability at the cooling front. This should be a testable prediction of this model, 
although this effect will be more dramatic in the later stages of the luminosity decline. 
Conversely, a signihcantly broader temperature transition prohle than the one seen here 
would be an indication of nonlocal ehects in a. 

In order to evaluate the size of the errors induced by using this solution, we need to 
examine the structure of the cooling wave near the cooling front. In this region the state 
variables will have a radial scale length much less than r, so the dynamical equations can 
be approximated by saving only radial derivatives of M and E. 

With this in mind we can write 


dtT.{r) ^ -rpdrT.. (34) 

Invoking equation (^) we have 

- rpdrT. = —drM. (35) 

27rr 

We can integrate this from E = 0 to S = E^ to obtain 

AM PS —ri7’27rri?Ei? = ^ Mp, (36) 

Vp 

where AM is the amount by which the mass outhow at the edge of the similarity solution, 
where E = 0, exceeds the mass outflow from the actual cooling front, where E = E^?. For 
the particular case examined by CCL we have 

AM 5.4 M^, (37) 

\rpiLpJ 

which will be somewhat less than half for (csp/rpQp) between 10“^ and 10“^. 

Although this amounts to a large correction to M, in agreement with the comparison 
between CCL and the prediction of the similarity solution, it does not imply that the 
behavior of M with time is subject to a similarly large correction. Assuming a Kramers law 
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opacity and n = 1.5 we can look at the evolution of the mass flow through the point where 
S = We hnd that 


tdt \n{MF - AM) 


tdtln 

= tdtMn - ( 


CsF 


-V 


rp'Q.p ) 

CsF 1 
\rVLp} 


= tdt In Mo 


1 - 




( 


cf \rF^F) 


(38) 


In other words, we get only a small error when we treat the time evolution of the mass flow 
through the S = Si? radius as though it were strictly proportional to the mass flow across 
the outer surface of the similarity solution. Since the realistic condition on the mass flux is 
that it is hxed at S = Sir, rather than at S = 0, the negative sign in the above equation 
actually implies that the realistic solution will give tdt In Mq slightly more than the one 
predicted from the similarity solution. For the calculation of CCL the discrepancy should 
be of order a few percent, which is not signihcant. 


Similarly we can estimate the error we make in equating the radius of the cooling front 
with the S = 0 surface by integrating equation (H) using equations (|^) and (|T^). We obtain 


Ar = rp 


3 CsF 
A rp^F 


(39) 


For the calculation of CCL this implies an error of a few percent. 


CCL pointed out that comparing the observed e-folding time for the soft X-ray 
luminosity in X-ray transients to the results of their simulation gives a dimensionless 
viscosity of 

CK ~ 50 



(40) 


This coefficient of 50 is surprisingly large, since one would expect that almost any theory 
which gave the correct scaling law would involve a coefficient of order unity. This estimate 
of the coefficient is roughly inversely proportional to the minimum hot state temperature, 
which is hxed by the physics of the opacity in the hot state. The only other obvious way 
to get a different result is by using a different value of A, the ratio of the outhow velocity 
at the cooling front to apCsF- However, we see from equation (^) that this enters into the 
expression for the speed of the cooling front only as ~ A^A_ Even if we replace the value 
calculated by CCL, A ~ 1/6, with A = 1 the estimate of the coefficient in the expression 
for a only drops to ~ 20. 


One last point we need to consider is under what circumstances we can ignore the 
structure of the cold region in discussing the velocity and structure of the cooling front. 
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From the equation of continuity we see that the column density in the cold disk region, just 
outside the cooling front, will be of order TjpiVp/ — rp). If this exceeds the maximum 
stable column density for the cold gas, then the solution described here is inapplicable, and 
the cold region will affect the velocity and structure of the cooling wave. The value of the 
maximum cold state column density is somewhat uncertain. Here we will use the results of 
Cannizzo (1993b). Assuming n ~ 3/2, we get 


'c,marr 


6.2 


Vo.iJ 


(41) 


where rio is the front radius in units of 10^° cm. We need this to be larger than Vp/ — rp 
which is 


-rp 


A 


50 

Olp ) 


(42) 


We see from these expressions that this condition is generally satished, although not by a 
large factor, with very little dependence on ai? or rio. Whether or not it will be satished for 
substantially different models for a or disk opacities depends on the specihcs of the models. 


In sum, the similarity solution described in this paper will give a good description of 
the progress of a cooling wave in a hot disk as long as four conditions are met. First, the 
outhow velocity across the cooling front needs to be substantially greater than the cooling 
front velocity. In practice, this means that the hot part of the disk needs to have a height 
to radius ratio of order 10“^ or less for the thermal instability associated with hydrogen 
ionization. If this condition is violated then the errors in the solution can become large. 
Other thermal instabilities will involve slightly different conditions, depending on the values 
of q and A. Second, the the opacity of the hot disk material, for a few e-foldings interior to 
the cooling front, has to be described by a power law. The explicit solution here assumed 
a Kramers law, but any power law approximation will yield a similar solution. Third, the 
a needs to be purely a function of local disk parameters. This last condition is signihcant, 
inasmuch as the only prediction for a which self-consistently allows for the effects of local 
magnetic instabilities and gives the correct scaling law is based on an explicitly non-local 
model (cf. [Vishniac, Jin, fc Diamond 19901 , [Vishniac fc Diamond 1992|) . We plan an 
extension of this work to cover this model and to see if an observably different cooling front 
structure can be recovered from it. Fourth, neglecting the structure of the disk outside the 
radius where rapid cooling sets in is appropriate only when the ratio of Vp to rp is not so 
large as to imply an unphysically large column density in the cold gas. This does not seem 
to be an important limitation for the systems considered to date, but this conclusion is 
sensitive to the details of the cold disk structure. 
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Fig. 1.— Surface temperature, column density, and mass flow as a function of radius for the 
cooling wave similarity solution using Kramers opacity and n = 3/2. The vertical axis is in 
arbitrary units, but M is normalized to 1 at small radii. 




